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1. Introduction 

Incoherent scatter radar (ISR) has become the most powerful means of studying the iono- 
sphere from the ground. Many of the ideas and methods underlying the troposphere and 
stratosphere (ST) radars have been taken over from ISR. Whereas the theory of refractive 
index fluctuations in the lower atmosphere, depending as it does on turbulence, is poorly 
understood, the theory of the refractivity fluctuations in the ionosphere, which depend on 
thermal fluctuations, is known in great detail. The underlying theory is one of the most 
successful theories in plasma physics, and allows for many detailed investigations of a num- 
ber of parameters such as electron density n 0 , electron temperature T e , ion temperature 
T; f electron mean velocity v e , ion mean velocity V{ as well as parameters pertaining to 
composition, neutral density and others. 

Here we shall review the fundamental processes involved in the scattering from a plasma 
undergoing thermal or near thermal fluctuations in density. We shall relate the fundamental 
scattering properties of the plasma to the physical parameters characterizing it from first 
principles. We shall not discuss the observation process itself, as the observational principles 
are quite similar whether they are applied to a neutral gas or a fluctuating plasma. These 
observational principles are dealt with in other ISAR presentations. 



Figure 1 . Volume element d(f) in volume V illuminated by plane wave. 


2. Volume Scattering from Random Irregularities. Continuum 

In the neutral atmosphere the scattering is derived from the following consideration: 

With an incoming wave: E 0 c~ % * in r [note that exp (hj 0 t) is understood] and a dielectric 
constant in medium: 


£ = c 0 + A‘(n), 


the polarization of the medium becomes: 


P( f)-Ae(r) 

This oscillating polarization acts as am equivalent current: 

j(r) = iuj 0 P(r) = iu 0 AzE 0 * ^ 


and the vector potential at a point far from volume V is: 


A{f) ~ 


4?r 







Observer 



Figure 2. Scattering geometry to calculate vector potential at the observer. 

with | r - r '[ % |r| - n\ ■ r\ where rii « 

We obtain: 

= [ A £ (r') e - r " ( *' i *-*"* , d(r') 

4tt |r] J v ' 

which means that the received field depends on the spatial Fourier component with 
vector k = kin - krec 
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Figure 3. Relation between the wave vectors of the transmitted and revised 
waves and the spatial Fourier component of the dielectric fluctuation. 

It is easy to show from this that for a plane incoming wave of flux (Poynting vector mag- 
nitude) Si n , the flux at the receiver is: 

Srtc = sin 2 X |Ae(fcj„ - k nc ) | a S.„ (6) 


ji 0 = 47T • 10" 7 Henry/m 
sin\ = polarization factor, see next section. 

Ri = distance between scattering volume and observer 


In a plasma 


A e(f) = - 


An(r ) • e 2 
m • cj2 


( 7 ) 


Where: 


An(r) = electron density fluctuation 

e = 1.602 * 10“ 19 Coulomb = elementary charge 
m = 9.110 * 10“ 3l J kg = electronic mass 

and one obtains by substitution: 

5 r «=(^-) \An{k in - k„ c )\ 2 ■ sin 2 x ■ S in (8) 

where r 0 is the classical radius of the electron defined by 
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0 4 tt me 2 * z Q 


6 o = 8.854-10 ~ u F/m 


c = 2.998 -10 s m/s 

3. Scattering from Individual Electrons. Discrete Particles 

Let us recapitulate the derivation of the scattering by a free electron of electromagnetic 
waves impinging on it. It was originally thought that the scattering from a plasma could be 
considered a super-position of scattering from individual free electrons and that the strength 
and spectral broadening can be used to determine density and electron temperature. For 
the ionospheric plasma and for the frequency used in such scatter experiments, it turned 
out that the scattering could not usually be considered that simply and that the actual 
situation was more favorable from an experimental point of view. Nevertheless, even with 
the more complex theory of interacting electrons the scattering from an individual free 
electron forms an important and essential ingredient. 

Assume as in Section 2 that the electric field set up at the position occupied by the electron 
is; 


f?(t) = £ 0 ■ e +lW (0) 

where E 0 is a complex electric field amplitude which allows for an arbitrary polarization. 
When we assume that uj 0 » where is the angular gyrofrequency of an electron, 
the equation of motion of the electron becomes: 

- E 0 e — m' v/e (10) 

Solving for u(<) with the substitution u(t) = v 0 one obtains: 

v 0 = -M — — E 0 (11) 

mw 0 

Note that we have neglected the spatial variation of the external electric field and the fact 
that the electron moves in this field. We have also ignored the force caused by the motion 
of the electron in the magnetic field of the incoming wave. Both of these effects could 
contribute to a ponderomotive force which we ignore. We also note that the motion of 
the electron is considered undamped. This cannot be strictly true since the electron, even 
without collisions, is re-radiating because of the oscillation and hence, must experience 
damping. 

The current density associated with the motion of this electron becomes; 


J(r, t) = -ev(t)6[r- r e (f)] 


( 12 ) 


where r t (t) is the position of the electron and where 6(f) is a spatial deltafunction. With 
such an oscillating current at the origin we obtain as in the previous section: 


A( r, t) = -«• £ 0 e +, <“'‘-*’ 

4irmuj 0 


r.' _ 

R i 


(13) 
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The position of the observer r may vary with time due to motion of the electron and we 
must substitute when necessary: 


r = r'(t) = 



where r l (t) is the radius vector describing the position of the electron relative to the position 
of the observer at time t. 


The Poynting flux vector at the receiver due to the radiation from this electron becomes: 

r. 1 | ff\2 1 ( e * \ l^rec x Eo\ 2 /i|\ 

S '“ = ’2 ,|fl| ~V\,l — (14) 

where 77 = yJii 0 /e~o = 376.7 ohms. Introducing the polarization angle \ through 


sinx 


__ l^rec X E 0 1 

and the pov^er density incident on the electron^ by: 

S in = \E 0 \ 2 h-\ 


(15) 


(16) 


we obtain: 


Srtc = (r a /Ri) 2 ■ ■ S in = 10 2 Stn 2 x • S in (17) 

The usual radar cross section a 9 is defined by: 

<7 S = 4 xR] • = 4 nrl sm 2 x — 10" 28 sm 2 x(ro 2 ) (18) 

*->«n 

In calculating the single electron scattering we have implicitly assumed that the driving 
electric field is linearly polarized and that the angle between the field and the direction of 
the receiver is x- I n actual fact the electron is often oscillating in two linearly polarized 
fields of arbitrary relative phase and amplitude. When this is the case the interpretation of 
sin 2 \ is more complex and the amount of scattered energy available may only be received 
provided the receiver is properly “matched” to the scattered wave. 

Consider a plane wave propagating along a positive z-axis. The complex amplitude E a may 
then be represented as; 

E 0 = Z o (cos0 ■ e z + e' 6 sinp ■ e y ) = £„ j j = £°P (19) 

A linear polarization along the x-axis corresponds to: 

Mi} 

and along the y-axis: 

Mi} 
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Circular polarization is represented by: 

right circular 
left circular 

It follows that the scattered field at the receiver can be expressed as follows: 

E r (r,t) = £-ni x (n, x (20) 

Hi 

where ni is a unit vector along the direction of the scattered wave. In order to specify 
the polarization both of the transmitted and scattered wave we introduce the following two 
coordinate systems: 

A. For the transmitted wave: 

€ z “ Tl 0 — 

e y = unit vector normal to plane defined by ^»nj ^rec 

e z = normal to both e y and e t in a right handed coordinate system. 

B. For the scattered wave: 

e z = Hi = ^rec/i^reci 



e' r = normal to e y and e f x in right handed coordinate system. 



»-U 


V5H' 



Figure 4, Relationship between polarization of incom- 
ing and scattered waves. 
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In terms of these two coordinate systems we may write: 


or: 



(21a) 


This determines the relationship of the two fields. In particular we note that 

sin 2 * = sin 2 P + cos 2 fi • cos 2 9 (22) 

where 9 is the angle between n 0 and n\. 

It is now clear that one may express the relationship between scattered and incident electric 
field amplitudes as: 


where: 


w = -jf^P 


* = 



p = 


{ ^ 1 } ( 8 eometry ) 
Zt- 

r cos0' \ . 

< ~ } (receiver) 

\ sinp'e' 6 ’ J V ' 


(216) 


4. The Scattering from a Collection of Electrons 

When many scattering electrons are present inside a volume V rather than a single one the 
observed field is given by the sum 

£r = ~^pY. e ' if,(t) (23) 

Rl 

where N is the total number of electrons within the volume V , where the polarization and 
geometry of all the electrons are the same so that we do not have to sum over different s 
The previous distance between a scattering electron and the observer, Ru now represents 
the distance from an origin within the scattering volume and r p (t) represents the position 
of electron number p within this volume. The vector k = k rec — k tn comes about as follows: 
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Figure 5. Scattering from electron at r p . 


In the scattering volume, apart from a phase factor, the incoming wave has the form: 

so that the field at electron number p is: 

= (24) 

From the previous section the field at the receiver due to electron p becomes 

E r = -- (o V - t~' kRlr (25) 

Ri 

If the distance from the receiver to the origin in the scattering volume is R\ we have: 

Rip & R\ ^rec * r p 

Substituting this into (25), ignoring an irrelevant phase factor (c lkRl ) and summing over 
all electrons gives eqn. (23) provided 

£=*r«c-*m, (26) 

We next imagine the particle density to be expanded in a spatial Fourier series through: 

"(r.O = (27) 

k 


where 


n(k, t) = j d(f)n(r, t)***' (28) 

Since the scattering electrons must be considered point particles at r p (t), p = 1,..N the 
density is: 
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N 


n(r, t ) = ^<5(r - r,(t)] 

P=1 


with a spatial spectrum 


N N" 

n(k, t) = X) e+<E ' <V( ‘ ) = E n *^’ *) (29) 

P»1 P=l 

A comparison with (23) shows that the observed scattered field from the electrons may be 
expressed as: 


S r = ~£- to *P ■ • n(ifc, t) (30) 

The complex amplitude of the received signal is 

Ar(t) = l ^°' M ' n( ^ <) (31) 
= A 0 n{k,t) 

< A * (t)Ar (i + r) > = A 2 0 < n* (k, t)n(k, t + r)> ^ 

= A 2 0 n 0 -V <n*(k,t)n r (k,t + r)> 

where < ... > denotes ensemble average, where individual particles are assumed to be 
independent, where n a is the mean electron density and where: 

< n* (k,t)n r (k,t + r) >=< >= p p (k,r) (33) 

is the autocorrelation of density fluctuations associated with a single electron. 

Similarly, the power spectrum received is determined from the Wiener-Khinchine’s theorem, 
and 


Prtc{u) - rto V ■ P at $ p (k,u) 


(34) 


where P ot is the power scattered by cm individual electron under the same geometrical and 
external conditions and where: 


*,(£«) = r°° < n; (j 4- r) > tT^dr = 

J — oo 


-i 

-L 


oo 

»+oo 


< e if[f,(l+r)-/,(l)l > . e -i<^J T 


oo 

+oo 


p p (k,T)e 


r dT 


We often need the integral: 

G(lc,u) = P p,(k, T)e~ iur dr 
With this function the spectrum of a particle becomes 


(35) 


(36) 





(37) 
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#,(f,w) = 2-Jk[G(*,w)]. 

It will become apparent in the next two sections that both the plasma response to an electric 
field as well as the thermal driving force can be described on the basis of the motion of 
non- interacting particles. It is, therefore, useful to study their motion in some detail. 

We assume that the static magnetic field B 0 is directed along the z-axis. It will be conve- 
nient to introduce “polarized” coordinates to describe the position of a particle: 

Cartesian: Polarized: 

r={x,y,z} r = {r +1 , r_!, r 0 } (38) 

The relationship between the two descriptions is 

r i = x = ^(r 1 +r_ 2 ) 

r -> - ds(*-'V) y = ^( r -i- r 2) 

t 0 = f z = r 0 (39) 

The advantage of these polarized coordinates becomes evident when we state the equation 
of free motion of a particle: 


(fog 

dt 


= —iaQv Q 



m 


(a = ± 1,0) 


(40) 


where q and m are the charge and the mass of the particle (electron or ion) respectively. 
From this we determine the relationship between the (past) velocity at time t 1 = t — t in 
terms of the present velocity (at t) 


v(t — r) = T(r) u(t) (41) 

where f(r) is a diagonal matrix with elements. 

[f(r)U = e’“ nr = ffa(r) (42) 

The past particle position can be determined similarly in terms of present position and 
present velocity through 

r(t — r) = f(t) — T(r ) u(f) (43) 

where T(r) is a diagonal matrix which determines the particle helical motion. The elements 
are given by 


iaQr _ i 

[T(t)U = ■ : qQ - =g a (r) (a = ±1,0) (44) 

The single particle autocorrelation, eqn. (33) now becomes: 

Pp (k, t) =< e i3e > (45) 

where 


<3or — ’ ff—a 


(46) 
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The actual form of p ? (k, r) depends on the statistical distribution of velocities, the magnetic 
field strength and the angle 7 between k and the magnetic field direction. 

For a Maxwellian velocity distribution: 

/ 0 (v) = ( 2*)-Wvr*e-W' (47) 

where 

v ih =T/m 

T = kinetic temperature (in energy units) 



__ y/2v t h 

Q 


(48) 


(49) 


It equals the ratio of the r.m.s. orbital velocity and the angular gyration frequency. 

We note that for weak magnetic fields and arbitrary 7, or for k || B 0 (i*e. 7 = 0) and 
arbitary magnetic field strength the autocorrelation becomes 

p p {k y r) (50) 

When k is exactly perpendicular to the magnetic field, the autocorrelation for the plasma- 
density becomes periodic with period T = 2?r/fi. The depth of modulation increases with 
R and decreases with the scale of the density fluctuation A = 2 t r/fc, see Figure 6. 
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For intermediate angles between k and B 0 the autocorrelation of the density fluctuation 
cam be regarded a product of two factors, see (48). 



-(**)* (?y«w a 7 

A) 

e V 2 / 

B) 

-(kR) 2 siTi 2 • sm J 


which are sketched in Figure 7. 



Figure 7* The two factors which in the autocorrelation of in the 
density fluctuation of non-interacting particles. 


The spectral function <t>(k,ui) becomes (see 37): 


<t> p (k,u>) = 2 Re |G p (£,u/)j 

(51) 

where 


G p & u)m J°° 

(52) 


The next section introduces particle interaction through a self-consistent electric field. 


5. Response to a Field Particle. Electrostatic I nteractions., 

We now introduce particle correlations in the plasma components (electrons, various types 
of ions) by consideri ng the particles smeared out into a continuum and by assessing the 
response of this continuum to each discrete particle separately. 

The plasma response to an electric field t)can be obtained by solving the Vlasov equa- 
tion to first order in 6n(r, v, t) and E(r y t). The perturbation solution for any of the species 
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is given by (see Appendix A) 


£n(r, t7,f) = - 


n 0 q 

m 



(53) 


Here E(r, t) is a small electric field which we shall later assume set up by the fluctuating 
charges, often referred to as the self- coftsi stent field. The integration is carried out along 
the unperturbed pair tide orbits which were studied previously. This is indicated by the 
primed quantities in the integrand. 

We make the following substitutions: 


t 1 = f-r 


- r)e-'^e+*- 


r(r)tf 


W r (A dfo(v') 

dv' ( 1 dv 


(54) 


Integrating over all possible v we obtain the induced density fluctuation: 

n(Jt, t) = M J~ dr E(k,t-r) j d{v)t (r) e‘ X 

We now take the Fourier transform with respect to time, and relate the electric field to the 
toted electric charge fluctuation Q(k^uj) by: 


(55) 


E(k,v) = ^jQ(k, u,) (56) 

We, therefore, only include longitudinal electrostatic interactions. This is exact when there 
is no external magnetic field. When such a field is present the longitudinal and transverse 
modes are coupled and the electrostatic approximation breaks down particularly at very 
long wavelengths. In diagnostic experiments on the ionosphere only short wavelengths are 
used and the electrostatic approximation is adequate. 

Substitution of (56) into (55) gives 

"(*,“) = <?(*>") J ~ dT ■ e_,wr / ‘W T >|? ei " < 57) 

For each of the species we now introduce 



Hence, we obtain for the fluctuation induced in species cr: 


n a (k,u) = — -x<,{k,uj) ■ Q{k,u) (59) 

We now have all the ingredients necessary to compute the density fluctuations in the plasma. 
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6. Calculation of the Fluctuation Spectrum 

Let us now assume that we have a plasma of electrons and one type of singly charged ions. 
For the electrons we substitute: For the ions we substitute: 


4-4 

2 2 

9* — -e 

q* -> 

me — i ► m 

m a — ► M 



a* a c 

a„ -* a.i 

foer * f o 

foe r -* Fo 

X& X* 

X«r X. 

n 9 — ► n 

n a - N 


Consider first the fluctuation associated with a particular electron . 

A. Intrinsic fluctuation n p (k, cj), See Section 4. 

B. Induced electron fluctuation nj(£, u) = 

(electron dressing on electron) 

C. Induced ion fluctuation N^(k,uj) = —X^Q t (lc,u)) 

(ion-dressing on electron) 


The total charge fluctuation associated with this single electron is 
Q c (k,u) = -e(n p + ^jQt) + e(-^-Q e ) 

from which: 


Q e {k,u) 


— e • n p (k, lj) 
1 + X« + Xi 


(60) 


(61) 


Since we are interested in the total electron-density fluctuation induced by electrons, not 
the charge fluctuation, we have: 


n e (k,u) = n p (k,w) + n\(k,u) 


= n P (k,u) + 


1 + X« + X. 


rcpC^X 1 + X.) 
1 + X« + Xi 


(62) 


The mean power spectrum associated with the thermal excitation by electrons is, therefore, 
found by averaging over the electron velocity distribution. If the velocity distribution is 
Maxwellian (47), then the independent electron spectrum is given by (51) and the result is 


. , ,r . < K(‘,“)l’ > 11 + Xil 1 _ 4.(^11 11 + ».P 

>= |i + *. + *,!■ U + X. + X.I 1 


(63) 


where $ e (fc, w) is the independent single electron mean power spectrum discussed in Section 
5. 
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Next consider the more important fluctuation arising from the thermal motion of an ion: 

A. Intrinsic fluctuation N p (k, uj), see Section 4. 

B. Induced electron fluctuation nj(fc,u;) = Qi(k,u>) 

(electron dressing on ion) 

C. Induced ion fluctuation Nl(k,uj) = — ^ Qi(k, u>) 


Solving for the charge fluctuation Qi(k,v) we obtain: 

eN,&u) 


Qi(k,U>) = 


1 + x« + Xi 

and the fluctuation induced in the electron density is given by: 


(64) 


nl(k,u:) 


u) 

1 + X« + X. 


n,(k,w) 


(65) 


The total electron density fluctuation in the plasma, therefore, can be expressed in terms 
of the independent particle spectra for electrons and ions and the response functions of the 
plasma as follows: 


< |n(£,w)| a > = [< |n e (£,w)| 2 > + < |n,(£,w)| 2 >\n 0 V 


( 66 ) 

H + y ; [ 2 $ t (fc,u>) + |y t l 2 $,(t,u>) 

|1+X«+X.| 2 

This is the form of the spectrum which has been widely used in the analysis of the incoherent 
scatter data, and which has remained valid for nearly 30 years! 

Fluctuations in ion density, charge density, electric field, currents, etc. can ail be obtained 
by an analogous procedure. An extension to a multi-ion plasma is relatively trivial. Differ- 
ent temperatures for electrons and ions are allowed. 

Collisions have not been considered but can, in some cases, be taken into account by 
regarding the particle motion as a stochastic rather than a deterministic process, see Section 
8. 


7. Discussion of Results 

The spectrum function contains the \t and together with and All of these 
functions are related to the autocorrelation of independent particle fluctuations as follows: 

Introduce for electrons (see 48): 




2 T' 
v ih = — 

m 

a -a = fc-o ' 9a 


D} = 


£ o T t 


{ga computed with Q = = 

(Debye length) 


m ; 



And for the ions: 


,„(*>) = Pl (fc,r) = e-W»MT 

T, 


v\ = — 
th M 

M — ionic mass 
A — or = k— a * <7 Or 

D] = 


eB a 


(g a computed with ft = H, = + -^-) 


2 


n 0 • e* 

With these definitions we obtain: 

2 


X«(^w) = ( — ) [1 +u> * /m jG e (fc, u;) j — iu? • i?e |G e (fc, u;)|] 
X,(jfc,w) = (1 + w/m |G,(fc,w)| -iw •iie|G i (fc,w)J] 


and: 


Here, from Section 4: 


$ e (J,w) = 2 -fle {<?,(*, w)} 

$,(fc,u;) = 2-i?e{G,(fc,w)} 

/ oo 

e - j 0 ?* ! J I J cos uiT dr 

/ oo 

with similar definitions for the ions. 

When there is no magnetic field or when k\\B 0 one obtains: 

ujRt {g(£,u>) j = Z e~> 2 

1 - z 2 + \Z K - JjZ* Z« 1 


1 + 


ulm |G(fc,u;)| = 


'ZT “ z* ~ ^ >> 1 


where we have put 
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Figure 8. Plasma dispersion functions. 

The only difference in the electronic and the ionic functions comes from the difference in 
thermal velocities. 


For values of u such that w < JcVtk x» and Xt are of the same order of magnitude. 
However, ~ RcG t j RtGi ~ V t h/vth ~ yJmjM 


It follows that as long as kD t and kD t << 1 the dominant excitation of density waves 
at low frequencies must stem from the ion excitation, see the numerator of (66). In the 
denominator x« — (iTJj) 2 whenever u> < kV i*. 

It follows that the low frequency part of the spectrum simplifies to: 


< Kt,u ,)| 2 


>« 2 n a V — i >— — 


(70) 


For T e = Ti the factor multiplying Re j<3i(u;,fc)J starts at J at u — 0. As T e /T, increases 
above unity the depression near w = 0 increases and a near line spectrum develops, see 
Figure 9. 

When kD t and kDi become much larger than unity, then the electronic part of the spectrum 
will dominate and we obtain a Gaussian spectrum with a width corresponding to the thermal 
motion of electrons. 

A resonance occurs near the electron plasma frequency w e . This can be established by 
looking for a zero in the denominator of (66) near the plasma frequency. Since x» has 
vanished near u> e we have: 


1 + 



{i 


+ uiIm(G t ) — iwi?e(G e )} 


= 0 


(71) 



350 



Figure 9 . The ionic part of the fluctuation spectrum for various values of 

As u ; is close to the expansion of the plasma dispersion functions for large 
cam be used and we obtain 



Where we have neglected the imaginary part The solution is the familiar expression: 

a ; 2 =wJ[l + 3(M) € ) a ] (73) 

the spectral peak associated with this oscillation is apparent in Figure 10* The peak can 
be strongly enhanced by the presence of photoelectrons or other suprathermal charged 
particles. The actual enhancement level involves the short-range Coulomb collisions as well 
as the angular distribution of the suprathermal electrons. 

Let us now briefly turn to the effect of the magnetic field. As far as the ions are concerned 
the gyrofrequency satifies the relation: 

kv th » n, 

This means that the correlation function 
becomes modulated with a periodicity TJ — 2n /H*. 

In practice, however, these modulations are blurred out because of the diffusion of the ions 
away from their deterministic orbits. Very close to perpendicularity with the magnetic 
field (i.e. 7 = 90°) the spectrum may become very narrow if the radius of gyration of the 


Te/Ti. 

arguments 
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Figure 10 . The fluctuation spectrum for T e /Tj = 1 for 
varying Debye length D e , X p = 

electrons is so small that the electrons are prevented from participating in the fluctuations 
of the ions. 

With a magnetic field the resonance associated with the plasma frequency becomes modi- 
fied. When >> Q ej when kR < 1 and when 7 is not too close to 90° one obtains: 

w 2 = w 2 [l + 3(fcD e ) 2 ] + ft 2 ■ sin 2 7 (74) 

An additional resonance, which cam also be observed, arises because of the presence of the 
magnetic field. As 7 — * 90° : 


Whenu;* > f one obtains: 




(75) 


u *i = oa (76) 

which is sometimes referred to as the lower hybrid frequency. The strength of these lines 
depend on the relative magnitudes of plasma frequency and gyrofrequency, on the presence 
of suprathermal electrons (or ions) etc. 

The total power residing in the electron plasma oscillation at thermal equilibrium is deter- 
mined by integrating the spectrum through the electron lines with the result that 


L 


$du> = 


electron line 


1 +(kD')* 


(77) 




1 


(78) 
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whereas the power in the ion line becomes: 


[ i + (kD t )*][l + lfi+(kD c y} 

Hence, whenever JcD t is small the ion contribution dominates. 

Note that (78) — ► \ when ( kD ) 2 — ► 0 or more generally —► y when T t ± T<. Hence, 

the “ion-scattering” is half of the free electron scattering, which is obtained from (77). 
When (kD) 2 — ♦ oo (77) — ♦ 1. 

8. The EfFects of Collisions 

Let us return to the solution of the inhomogeneous Vlasov first order equation: 

6n(r,v,t) = -^ (79) 

Remember that the past R, V at time t - r approach r, u at time f, and that the posi- 
tion/velocity travel back in time exactly as in deterministic orbit in accordance with the 
equations of motion of free particles. 

Suppose that the particles suffer collisions. In this case, it is no longer possible to say for 
certain where they came from, because the previous history of the particle arriving at r, v 
at time t is a random process. If the collisions occur with like particles the process becomes 
difficult because pairs of particles have unrelated motions in this theory. However it is 
often the case that the collisions occur with particles of another kind with little dynamical 
mutual coupling. Important examples are: 

ions in low ionosphere collide with neutral gas molecules which 
are much more numerous than the ions. 

electrons deviating from their deterministic orbits because they 
have to move in the random field of near-stationary ions. 

In cases such as the two quoted we introduce: 

W+(r, v, t\f - R, \f y t - r) = conditional probability density 

of finding a particle at r, v 

at time t given that at time t — r it was 

at r — R,V 

The joint probability density of the two events (r, u, t) and (f— R, V, t — r) is 

W + ( r,v,t\r-R,V,t-r)f 0 (?) 

The individual probability of r - R, v, t - r given that the present coordinates are r, v,t is: 

W_(r — R,V ,t — r|r, u, t) 

and the joint probability density of r, v, t and r — H, V, t — r is: 




353 

W.(f-R,V,t-r\r,v,t)f 0 (v) = 

= W+ft v, t\f- R,V,t- r)f 0 (9) (80) 

but, 

W-(-R,V y -T\o y v,o) 

clearly must equal 

W+(R, V,t\o,v,o) 

from the symmetry of the equations of motion. 

Taking the spatial Fourier transform of the perturbation solution one obtains: 

n(£,v,t) = -— re ( t,t-T)<mp .'***> dr (81) 

m J 0 ov * 

where the average is taken over all the different particle orbits which lead to (r, v, t). 
Explicitly with ■ f 0 : 

n(f, v,t) = + -^ J~ i(l it- t)JJ W + (R t V,r\o,v,o)V ■ e +itR d(A)d(9) ■ f.{v) (82) 
Hence, if we introduce: 

fo(v) jj d(R)d(V) W + (r)V ■ e+ iiR = A(k, v, r) (83) 

we obtain: 

n(*. «?, <) - -gj J r £(£, t — r) ■ A(k, v, r)dr (84) 

Let us formulate the modified scatter theory in terms of transition probability averaging. 
Introducing as before: 

i(?, w ) = ^gM (ss) 

and taking Fourier transforms, one obtains: 

= f dTe ~ iwrJ ' (86) 

Integrating over all arrival velocities one finally obtains: 

n(f,w)---x(E.«*0 <?(*,«) (87) 

Q 


where now 
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x(k, U,) = -i jf * e "' wr / «*) • /.(«> 

JJ d(&)d(i 0 W%(£v,r|o,tT,o)(BV ,J * 

Similarly - studying independent particles - diffusing along as a result of collisions 
another gas - but ri£t exposed to an external field - one obtains the expression: 

\n p (k,u\ 2 =2 J° dr • e~' UT J d{v) f a {v)- 

J I d(R)d(v ) W+(R, V, t\o, v, o) t~' iR 
So that what we previously referred to as G now takes the form: 

G(£ t u) = jf° dr t~ iur J d(v)f 0 (v) JJ d(R)d(V)- 

■W + (R,V,r\o,v,o)e- iin 

By properly manipulating the expression for x(^ w ) we obtain 

1 




■ ( 1 + c o IttlG — iuJ Re G) 


(kDy 

which is of the same form as before. Consider a model for W(R y v, r|o, v,o) : 
Suppose the particle is moving as if in Brownian motion with an equation 

— = —fiv + A(t) (Langevin's equation) 
dt 

(do not confuse with 0 in Section 2!) 

Then, from Chandrasekhar’s work: 

W + (R,V,t\v) = 


8 ^(FG - H*) 3 / 2 


a R*-ihA c 5+f S* 
e ~ 7{ra-B*) 


R 0 = R-§(l-e-» r ) 

S = V-v-e~ fir 

F =^[2^r-3 + 4 -e-^-e- 2 ^] 

G-St 1-^1 

The correspondence for /? — * 0 

W+(R,V,t \v) = 6(R-vt) 6(V -v) 


(88) 

with 

(89) 

(90) 

(91) 

(92) 

(93) 


(94) 
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This is a formulation which preserves particles, momentum and energy (with the back- 
ground). We shall now use it for a study of plasma line enhancement due to photoelectrons, 
or other energetic “tail” electrons, neglecting the magnetic field. 


By substitution one obtains: 

. e - iu *dT 

If we assume 0 to be small we obtain: 



G(k 


',«) * j r t- ar7 - iuT dr + ^ jf° r 3 e~ aT * ~' UT dr 


Go(k^) 


with a = — 


(95) 


(96) 


Hence: 


G(k,uj) = G a (k,u>) - i -j G 0 (k,u) 


^/friction term collision frequency 

(97) 


The asymptotic expansion for Im G for large — Z is 

-- {1 + Z~* +3Z~*...} 

w v 

Hence, we obtain approximately: 

G(k, w )«G 8 (f, w ) + ^-2-3-3j... 


k 2 v\ 


th 


+ 


(98) 


= G 0 (k,u) + P- 

=c „ (J , u)+ ,.% i 

Turning to the fluctuation spectrum for electron plasma oscillations one obtains 


(99) 


< Hk,uj)\ >- 


u +x«i 


(100) 


l + Xe = 1 + 



+ u) Im G — uilm G — iuiRe G] 


w 1 £ — 3 


(kvo ,) 2 




(kDy 


u) ■ Re G 


Denominator: 
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» i - 4 - 3 (^) 2 ■ 4 


Expand to first order in u about w = u? r 

1 + x '* + ^ (w ~ Wr) “(W w * eG 

= _( w _ Wr ) + __.v/^— (e i z ) + ... 

The power spectral shape in this approximation becomes: 

2n 0 Re G 


<\n(k, w)| 2 >~ 


^(w-w r ) J +(iD) 4 (37-) J 


( 101 ) 


( 102 ) 


(103) 


This is a Lorentzian which can be integrated. If this is done one obtains for the intensity: 


1 fmjv£) + /,.(*>) ±£ 

WDy /-(»♦) + 


(104) 


When electron-ion collisions dominate, /? can be expressed approximation by as u p 
(A ~ number of electrons in Debye sphere.) 




Wj 

* 


fm ~ n m ( 2 ^)^ 2a ’ t~ mv +f 2T background plasma 
f p ~ same for hot plasma. 

Collisions between ions and the neutrals cause the mobility of the ions to decrease. The 
collisions, therefore, effectively damps the ion-acoustic waves, which causes the frequency 
spectrum to narrow. Studies of the narrowing of the frequency spectrum with decreasing 
height can be used for studies of the neutral density at D and the ^-region heights. 

9. S u m mary apd some final remark ? 

A summary of the application of the ion line of the incoherent scatter radar technique to 
ionospheric measurements can best be given in terms of the sketches shown in Figure 1 1 . 

The plasma line, in addition to providing information on suprathermal flux, see equation 
(104), also provides the possibility to determine the electron temerature and the electron 
density. Ignoring the geomagnetic field for simplicity we have, for the effective vectors k x 
and i 2 (= — k rec ) the following plasma line frequency offsets apply, see equation (73): 


/i 2 fl = /« 2 (l + 3(fc 1 i?) 2 ] 


/ 2 2 k = /« 2 [1 + 3(M>) 2 ] 
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Figure 11. Sketches showing various observations leading to several plasma parameters. 


If the two different fc- vectors are generated by shifting the radar center frequency, the 
plasma line spectra at the two radar frequencies /oi and /oi will appear as shown in Figure 
12. Different Jfc-values can also be obtained by changing the geometry in a bistatic setup. 
In the case shown in Figure 12 we have: 
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(/ 2 ft - /.ft) (/2ft + h ft) = 3Z) 2 • f\{k 2 - k x ){k 2 + k x ) (107) 

-V * 

~ 2 /. 


, , * wi ft 47t(/ 02 - /oi)4tt(/o 2 + /oi) 

0/ = hR- }\R — o D y — ^2 " 

(108) 



Figure 12. Plasma line spectra at radar frequencies /oi and /o 2 - 

All of the quantities involved are well determined except for T t which can be found this 
way. 

Numerical example (EISCAT parameters) 

T e ~ 2000 ° K 
f 02 = 933 MHz 

foi = 224MH z 6f % 335 kHz 

f t = 5 MHz 

Hence, from a two-frequency plasma line experiment one can deduce the electron temper- 
ature accurately and independently. 

In Section 2 of these lecture notes we started out considering the scattering from random 
irregularities in the dielectric function. Throughout the dielectric considered in this pa- 
per was a nearly lossless plasma. However, a neutral gas also exhibits random density 
fluctuations, and the curious readar might wonder whether they are detectable at radio 
wavelengths. In order to answer this question, consider a gas dense enough to support 
sound waves. The density fluctuation may, therefore, be considered as a superposition of 
thermally excited sound- works of varying wavelength and direction. As in Section 4 we ex- 
pand the parameters (velocity, density, temperature) associated with the acoustic wave-field 
: n a spatial Fourier series: 
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rl ( r " < ) = 7? n( *’ <) ' e ”’ X ' 

k 

where n can be density, pressure, velocity etc. The wave-field amplitudes n(fc, t) all must 
satisfy the wave-equation 

n(k,t) + k 2 c]n(k,t) = 0 (109) 


where 



( 110 ) 


where T is the mean gas temperature in energy units, M the molecular mass and 7 (1.4 for 
air) the ratio of specific heat at constant pressure and constant volume. From equipartition 
arguments (assuming minute losses to insure equipartition) as used in deriving specific heat 
of solids we find that 


< KM p >= 


[<5(u> — keg) -f 6(u> 4- &c,)] 


( 111 ) 


FYom Torn Sato’s lecture, his equation 2, we see that for dry, nonionized ain 

1.55 • 10~ 9 v(k)(mb) _ 


A e(k) = £ 0 


T(°K) 


= £0 


1.55 • 10~ 7 v(k)(N/m 2 ) 
T(°K) 


(112) 


= £ 0 * 1.55 * 10~ 7 £ * k * n(fc) 

where 

k — Boltzmann's constant = 1.38 * 10“ 23 J/° K 
combining 111 and 112, one obtains: 

< |A£(Jt)| 2 >= £jj 7 • k 2 • n 0 • V ■ 2.4 • 10“ 14 (113) 

The radar cross section per unit volume of the gas is formed by combining equations (6) 
and (18) to give: 

<r = ^2. 3.73 . io-«° n 0 (m 2 /m 3 ) 

4 IT 

As a specific numerical example, take the Arecibo S-band wavelength and the atmospheric 
number density at sea level* 


A 0 = 0.125 m 


n 0 = 2.688 *10 28 m“ 3 

This gives for the troposphere: 


vtrof « 5 • 10 26 m 2 /m 3 
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As a comparison the cross section per unit volume in the F-region is typically: 

<7 F « 10" 18 m 2 /m 3 

If the Arecibo beam could be focused at lower altitude, if full advantage could be taken of 
the reduced distance and the reduced bandwidth of scattering one might be able to make 
up for the nearly seven orders of magnitude discrepancy in specific radar cross section. It 
is probably a much more practical approach to excite low frequency sound waves as done 
with the Shigaraki MV radar system, and scatter from them. 



